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Abstract 

We provide a detailed treatment of relativistic Lotka-Volterra hierarchy and a 
kind of initial value problem with special emphasis on its the theta function rep- 
resentation of all algebro-geometric solutions. The basic tools involve hyperellip- 
tic curve K n associated with the Burchnall-Chaundy polynomial, Dubrovin-type 
equations for auxiliary divisors and associated trace formulas. With the help of 
a foundamental meromorphic function <fi on K, p and trace formulas, the complex- 
valued algebro-geometric solutions of of RLV hierarchy are derived. 



1 Introduction 

Nonlinear integrable lattice systems have been studied extensively in relation 
with various aspects and they usually possess rich mathematical structure such as 
Lax pairs, Hamilton structure, conservation law, etc. The Toda lattice (TL), 

(1.1) 




is one of the most important integrable systems. It is well-known soliton equations 
such as the KdV, modified KdV, and nonlinear Schrodinger equations are closely 
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related to or derived from the Toda equation by suitable limiting procedures (61IH1- 
Another celebrated integrable lattice system is Lotka-Volterra (LV) lattice, 



LV : 



Ut = u[v — V 
v t = v{u + — 



u), 



(1.2) 



or 



a t = a{a 



+ 



a 



(1.3) 



by setting a(2n — 1) = u(n), a(2n) = v(n). Ruijsenaars found a relativistic inte- 
grable generalization of non-relativistic Toda lattice through solving a relativistic 
version of the Calogero-Moser system d.The Lax representation, inverse scat- 
tering problem of the Ruijsenaars-Toda lattice and its connection with soliton dy- 
namics were investigated. A general approach to constructing relativistic general- 
izations of integrable lattice systems, applicable to the whole lattice KP hierarchy, 
was proposed by Gibbons and Kupershmidt @. However, nobody have known 
what the relativistic Lotka-Volterra lattice is until Y. B. Suris and O. Ragnisco 
found itfll. 

There are very close connnection among the Toda lattice, the Lotka-Volterra 
lattice, the relativistic Toda lattice and the relativistic Lotka-Volterra lattice. It is 
well known that ( 11.11) and ( 11.21) are both discrete version of the KdV equation in 
the sense of different limiting process and in the non-relativistic limit h — > 0, the 
relativistic Toda system (RT), 



and the "relativistic splitting" relativistic Lotka-Volterra systems (RLV), 



reduced to the well-known Toda lattice equation (11.11) and Lotka-Volterra lattice 
equation (11.21) , respectively. Moreover, the Miura relation can be summarized in 
the following diagram: 




(1.4) 




u t = u{y — v + huv — hu v ), 
v t = v{u + — u + hu + v + — huv). 



(1-5) 



RLV 



LV 



RT 



TL 



2 



RLV 



n 



RT, RLV 2 



RT 



Here 




a = uv 



b = u + v 




a = u^v, 
b = u + v. 



Mathematical structure related to relativistic volterra lattice (11.51) such as Lax 
integrablity 0,2x2 Lax representation ifTO) . conservation laws ifTTI . bilinear 
structure and determinant solution lfT2) have been closely studied and hence the 
purpose of this paper is to uniformly construct algebro-geometric solutions of the 
relativistic Lotka- Volterra hierarchy which invariably is connected with geometry 
and Riemann theta functions, parameterized by some Riemann surface. Algebro- 
geometric solutions (finite-gap solutions or quasi-period solutions), as an impor- 
tant character of integrable system, is a kind of explicit solutions closely related to 
the inverse spectral theory |[T5l,[T7) . Around 1975, several independent groups in 
UUSR and USA, namely, Novikov, Dubrovin and Krichever in Moscow, Matveev 
and Its in Leningrad, Lax, McKean, van Moerbeke and M. Kac in New York, and 
Marchenko, Kotlyarov and Kozel in Kharkov, developed the so-called finite finite- 
gap theory of nonlinear KdV equation based on the works of Drach, Burchnall and 
Chaunchy, and Baker |fT3l[T6*) . The algebro-geometric method they established al- 
lowed us to find an important class of exact solutions to the soliton equations. As 
a degenerated case of this solutions, the multisoliton solutions and elliptic func- 
tions may be obtained ifTTH Tl. Its and Matveev first derived explicit expression of 
the quasi-period solution of KdV equation in 1975 Ifl4) . which is closely related 
to the finite-gap spectrum of the associated differential operator. Further excit- 
ing results appeared later, including the finite-gap solutions of Toda lattice, the 
Kadomtsev-Petviashvili equation and others (51 [TH |7), which could be found in 
the wonderful work of Belokolos, et al lfT7) . In recent years, a systematic approach 
based on the nonlinearization technique of Lax pairs or the restricted flow tech- 
nique to derive the algebro-geometric solutions of (1+1)- and (2+l)-dimensional 
soliton equations has been obtained ||25l |26) . An alternate systematic approach 
proposed by Gesztesy and Holden can be used to construct algebro-geometric so- 
lutions has been extended to the whole (1+1) dimensional continuous and discrete 
hierarchy models JIl[IS|27l|28l|29). 

The outline of our present paper is as follow. In section 2, the relativistic 
Lotka- Volterra equation (11.51) is extended to a whole hierarchy through the poly- 
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nomial recursive relation. The hyperelliptic curve associated with RLV hierarchy 
is given in terms of the polynomial. In section 3, We focus on the stationary RLV 
hierarchy, based on the polynomial recursion formalism introduced in Section 
2 and a fundamental meromorphic function <fi on the hyperelliptic curve K n , we 
study the Baker- Akhiezer function trace formulas, from which the algebro- 
geometric solutions for stationary RLV hierarchy are constructed in terms of Rie- 
mann theta functions. In section 4, we extend the algebro-geometric analysis of 
Section 3 to the time-dependent RLV hierarchy based on a kind of special initial 
value problem. Finally, in section 5 we give Lagrange interpolation representation 
that will be used in this paper. 

2 Relativistic Lotka-Volterra hierarchy and associ- 
ated hyperelliptic curve 

In this section, we investigate the relativistic Lotka-Volterra hierarchy and 
then derive the hyperelliptic algebraic curves associated with the algebro-geometric 
solutions of the newly constructed hierarchy. Throughout this paper, we have the 
following definition. 

Definition 1. We denote by £(Z) the set of all the complex-valued sequences 
{f{ n )}n=~oo- This is a vector space with respect to the naturally defined opera- 
tion. A subspace £ 2 (Z) C £{Z) is defined by the set of {f G £{Z)\ j2n=-oo \f( n )\ 2 < 
+oo, n G Z}. 

Definition 2. We denote by S* 1 * 1 the shift operators acting on ip = {ip^n)}^^^ G 
£{Z) according to (S ± ijj)(n) = ip{n ± 1). The identity operator I acting on 
if) = {^{n)}n=-oo G £{1) according to (J^)(n) = ip{n). We also define ^ = 
S^ij) G £{Z). 

We introduce the following 2x2 matrix problem |[T0l 



where ^ = (#i,^ 2 ) T G ^(Z) x £(Z),f = (p-l)^ 1 ^ = p,q = e (S+-r) 
and p, q, p, q axe potential functions and A is the spectral parameter. Here p = 




(2.1) 



4 



p(n,t),q = q(n,t) G £(Z),(n,t) G Z x R in time-dependent case and p = 
p(n), q = q(n) G l(Z),n G Z in stationary case. 

Define sequences {af(n)}* 6 N , {h(n)}ien and {cf(n)}* 6 N in£(Z) recursively 

by 

pa e - pa t = a e+1 - a e+1 + fb t - qc e , I G N , (2.2) 
- ~ <7 a 7+i = <ia>£+i + h, t G N , (2.3) 
fa7 +1 + c7 = pc e - Q + i - fa/+i, i G No, (2.4) 
a = 1/2, 6 = ~<? + , c = -f. (2.5) 

Explicitly, one obtains 



ai = -g + r + Si/2, 

bi = —p + q + + (g + ) 2 f + g + g ++ f + + q ++ — q + Si, 



C\ = —pr + q + r 2 + f + qrf — r8\ 



(2 6) 

a,2 = —p + q + f + (g + ) 2 f 2 — pg + f + g ++ f + q + f~ + g ++ g + f + r + q + qrr~ 
- g + r5i + «y 2 /2, 



Here {Si}g e ^ denote summation constants which naturally arise when solving 
(l2T2l)-(l231). 

If we denote by d e = a e \ S:j= o tj=1 ... t t,b e = b £ \ Sj=Qd=1 ... ti ,c i = c £ \ s . =od= i r .j, 
£ G N the homogeneous coefficients of b e , q, that is, 



cl\ = —qr, 

bi = —p + q + + (q + ) 2 r + q + q ++ f + + q ++ , 

C\ = —pr + q + f 2 + f~ + gff~, (2.7) 
a2 = —p + q + f + (g + ) 2 f 2 — pg + f + g ++ f + q + f~ + g ++ g+f + f + q + qff~ 



By induction one infers that 



5^_ fc a fc , £gN , (2.8) 



fc=0 



h = J2 6 ^h, ^N , (2.9) 



fc=0 
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i 

ct = ^2Si- k Ck, £eN , (2.10) 



fc=0 

introducing 5o — 1. 

To construct the relativistic Lotka-Volterra hierarchy we will consider the fol- 
lowing ansatz 



V (A) = ( J ^ 2n + 2 ^ ^2n+l(A) 
V <^2n+l(A) -D211+2W 



n E N , (2.11) 



where A 2n+2 , -B 2n+1 , C 2ri+ i and D 2n+ 2 are chosen as polynomials, namely 

n+l 



A 2n+2 (X) = aA~ {2n+2 - 2e) = a Q \-^ + ai \- 2n +■■■ + a n+1 , (2.12) 

n 

B 2n+1 (X) = b £ \- {2n+1 - 2e) = M~ (2 " +1) + h\-^ +■■■ + b n X-\ (2.13) 
1=0 

n 

C 2n+1 (X) = Y,^- {2n+1 ~ n) = c A- (2n+1) + c x X~^ +■■■ + c n X~\ (2.14) 



n+l 



D 2n+2 {X) = ^^A- (2re+2 - 2 ^ = rf A^ 2 " +2 ) + rfxA" 2 " +• ■ ■ + d n+1 (2.15) 



and the coefficient {a^}"^ 1 , {&£}" =0 , {q}"=o> {^}"=o ^ defined in £{Z). 
The stationary zero-curvature equation 

= U(X)V n (X) - V+(X)U(X) (2.16) 

is equivalent to the following equalities 

(pX - A" 1 )^, + qC 2n+1 = (pX - A -1 )A 2n+2 + rB 2n+1 , (2.17) 

(pX - X~ l )B 2n+l - qD 2n+2 = qA 2n+2 + 5 2re +iA, (2.18) 

?A 2n+2 + C 2 - +1 A = (pA - X~ l )C 2n+l - rD 2n+2 , (2.19) 

fB 2n+1 - D 2n+2 X = qC 2n+ i - D 2n+2 X. (2.20) 
From (12.161) , one finds the matrix V n {X) is similar to V+(X) and then we have 

trace (K(A))= trace {V n + (X)) , 
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namely, 

^2n+2 — -^2n+2 = ^2n+2 ~ -D2n+2- 

Hence A 2n +2 — D 2n +2 is n-independence. Without loss of generality we can 
choose 

d £ = a £ , I = 0, 1,2,- • • ,n, d n+1 = A, (2.21) 

where A = 1 = (■ • • , 1, 1, 1,- • •) e £(Z) is a constant value sequence. In station- 
ary case we have 

A.2n+2 = D 2n+2 , (2.22) 

since this can be always be achieved by adding a constant coefficient polynomial 
times the identity matrix to V n , which not affect the stationary zero-curvature 
equation (I2.16I ). 

Plugging the ansatz (12.12M2.15b (12.21b into (12.17b - (12.20b and comparing the 
coefficients yields the following relations for {a^}"= , {be}g =0 , {ci}™ =0 . That is, 



pa e - pa e = a e+1 - ae+i + rbe - qc e , I — 0, 1,- • ■ , n - 1, (2.23) 

pb~ t - b~ i+1 - qaj +1 = qa e+1 + be, £ = 0, 1,- • • , n - 1, (2.24) 

raj +1 + cj = pa> - ct+i - rae+i, I = 0, 1,- • ■ ,n - 1, (2.25) 

f&7 +1 - aj +1 = qc t+1 - a i+1 , £ = 0,---,n-l, (2.26) 

pa~ +1 - pa n+1 = 0, (2.27) 

pb~ - q A = ga n+ i + b n , (2.28) 

ra~ +1 + c~ = pc n - fA, (2.29) 

- A" = -A. (2.30) 



Hence if {at)e=o, {^}™ =0 , {q}£ =0 are defined by (|2T2l -(l2"3T). then one finds (12^231) - 
(12T25T) (12301) hold naturally. From (l2!23l - (|2!25l) . one calculate 

rb e+1 - qa e+ i 

= r (pbj - qaj +1 - qa e+1 - be) - q (pc e - ra e+1 - fa i+1 - cj) 

= p(rbj - qc £ ) - rb e + qcj (2.31) 

= p(aj - at) - rbe + qcj 

= a e+1 — O't+ii 

where we used the induction 

rb k ~ = qc k - a k , k = 0, !,■■■,£ 
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in the third equality of (12.31b . Therefore (12.26b is the direct result of (12.23b - (12.25b . 
Using (I2~2b-(l231). (l2T27T )- (l230l) are equivalent to 

pa~ +1 - pa n+1 = 0, (2.32) 

K+i + 9 a n+i - gA = 0, (2.33) 
c n+1 + ra n+1 + f A = 0. (2.34) 

Noticing the condition (12.26b and above analysis, ( 12.32b and (12.33b give rise to the 
stationary relativistic Lotka-Volterra hierarchy, which we introduce as follows 

s-RLV n (p, q) = s-RLvC(p, q, r) = ( a <£) L " ° n+1 ) = 0, p G N . 

(2.35) 

Explicitly, 

s-RLV (p, g) =s-RLV (P,?,r) 

— gf~ + q + f \ 
-pq + qq + f + q + -q(A + 6 1 /2) J (2.36) 

pq-p~q~ + q~ -9 \ n 

where we use the relation 

q+/q = e (S+-ir^ /e (S+-l)-\ = q ^ 

qr = p — 1 = p — 1, (2.37) 
gf + = (p — l)g. 

Taking <5i = — 2A = —2, (12.36b gives rise to the stationary relativistic Lotka- 
Volterra equation 

/ pq — p~q~ + q~ — q 
pq — p 



0, (2.38) 



which means 



pq — p q + q — q 
p+q+ — pq — p+ -)_ p 

I \ / pq — p~ ~q~ + q~ — q \ (2.39) 
S + — I J \ pq — p 

0. 
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Next we turn to the time-dependent relativistic Lotka-Volterra hierarchy. The 
zero-curvature equation is 

U tn (X) + U(X)V n (X) - V+(X)U(X) = 0. (2.40) 

C/(A),K(A),{a,(-,t)}^ by 
(ITTi d2j)-d22]), (I2TTT1) - (|215T) and (127ITT) . Then (l274"0l) implies 

{~Z X TH£ £)=<>■ 

where 

A n = (pX - X^)A M + qC 2n+1 - (pX - X~ 1 )A 2n+2 + rB 2n+1 , (2.42) 



A12 = ipX - X~ [ )B 2n+1 - qD 2n+2 - qA 2n+2 - B 2n+1 X, (2.43) 

A 21 = rA~ n+2 + C 2n+1 X - (fiX - \~ 1 )C 2n +i + fD 2n+2 (2.44) 

A 22 = rB~ n+2 - D 2n+2 X - qC 2n+1 + D 2n+2 X. (2.45) 
Inserting (l272l)-(l23T) into (12.411) then yields 

- Pt n + P(a~ +1 - a n+l ) = 0, (2.46) 

- Qt n + K+i + WZ+i - ?A = (2.47) 

- f tn + c n+1 + fa n+1 + f A = 0. (2.48) 



Using (ESI) (I2T261) (12461) (12747b . we directly calculate 

h n = d tn {{p- i)/q) =PtJq- (p- i)<k/<f 

= P(«n+1 - On+l)/q- (?» - 1)(K+1 + Wn+l - <?A)/<f 

= P( a n+1 - a n+l)/q- r{PK ~ b n~ ?«n+l ~ gA)/<? 

= p( a n+i - a n+1 ) / q - fpb~ / q + fb n / q + (fa n+1 + fA) 
= [p(a~ +1 -a n+1 ) + f(6 n (f<Wi + fA) 

= \P( a n+l ~ a n+l) + f(~K+l - Wn+l ~ <3<+l]/<? + (™n+l + f A) 

= - + (1 - - On+l) ~ fb~ +1 ]/q + (™n+l + ^ A ) 

= K+i - a n+i - r&~ +1 ]/g + (fa n+1 + fA) 
= c n+ i + fa n+ i + fA. 

(2.49) 
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Thus, ( |2.46l )- (|2.48l ) are essentially two independent equations ( 12.461 ) and ( 12.471) . 
Varying n E N , (12.461) and (12.471) give rise to time-dependent relativistic Lotka- 
Volterra hierarchy, which we introduce as follows 

RLV n (p, q) = RL\C(p, q, f) = ( f\ +/ (fl ^ ~ an+l ], ) = 0, p G N . 

(2.50) 

Explicitly, 

RLVo(p,g)=RLVo(p,g,r) 

-Pto +P(-qr~ + g + f) 
-9to + + q + - g(A + <Ji/2) 

-Pto + P(P9 - + <T - 9) 
^ (e( 5+ - 7 )" lln9 ) + e ( 5+ - 7 )" lln 9 (pq-p- (Sj/2 + A)) 

-p*o + p(p? - p~q~ + q~ -q) 

e (S+-/)- 1 ln 9 ( 5 -i _ I)~\ qt J q ) + e (S+-/)- 1 in 9 { pq -p- fa/2 + A)) 
= 

(2.51) 

where we used 

^ (e (S+-/)-ln 9) = e (^-/)- lng ^ ((5+ _ /) -l lng) = e (^)-^ q{S -l_ ir l {qt 

(2.52) 

in the third equality of ( 12.5 lb . Taking 5\ = — 2 A = —2, (12.511 ) is just the relativis- 
tic Lotka-Volterra equation, which is equivalent to the form of 



RLV (p,g) =RLV (M,f) 

-pto + p(p? - p~<T + <T - g) 
-?*o + qip + q + -p + -pq + p) 



(2.53) 



0. 



In order to derive the algebraic curve associated with the relativistic Lotka-Volterra 
hierarchy, we need the some change about U(X), which is 

U(\)^U(Z) = ZU(C 1 )=( P ~£ Z f ), * = (2.54) 

Let 

V n (0 = V^r 1 ), (2-55) 
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then U(£), V n (£) remain satisfy the zero-curvature equation (12.161) and ( 12.401 ), that 
is 

U(OV n (0 - V+(Z)U(0 = (2.56) 

K(0 + U(OV n (0 - V n + (OU(0 = (2.57) 
and hence the final form of (12.351) (12.501) is invariant. Assume 

*® = = -fil)) <2 - 58) 

and then one derives 

A n+1 (z) =A 2n+2 (C 1 ), (2.59) 

R n (z)=C 1 B2n+i(C 1 ), (2-60) 

C n (z)=r 1 C 2n+ i(r 1 ), (2-61) 

D n+1 (z) = D 2n+2 (C l ) (2.62) 
by comparing with (12.1 II ). Then (|2.17t - (|2.10t change into 

(p - + z<?(7~ = (p - 2)i n +i + z?B n+ i, (2.63) 

(p - z)fl" - qD~ +1 = qA n+1 + B n , (2.64) 

fi~ +1 + (7" = (p - z)(7 n - rD n+1 , (2.65) 

zf£~ - = zqC n - D n+ i, (2.66) 

respectively. 

Taking into account (12.561) (12.221) . one infers that the expression R 2n+2 (z), 
defined as 

R2n + 2(Z) = -(in+l) 2 ~ zB n C n , (2.67) 

is a lattice constant, that is, R 2n+2 = R 2n+2 , since taking determinants in the 
stationary zero-curvature equation (12.561) immediately yields 

p(l - z)[-(A n+1 ) 2 - B-C- + (A n+1 ) 2 + B n C n ] = 0. (2.68) 

Hence, R 2n +2 only depends on z, and one may write R 2n +2 as 

2n+l 

R2n+2(z) = -(1/4) H(z- E m ), n e N . (2.69) 

m=0 
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Remark 3. (i) Taking the transformation p — > hp, q — >■ hq, one concludes that 
the equation t\2.35\) and l\2.50\) change into the normal form of relativistic Lotka- 
Volterra system ITOl/ . 

(ii) Taking the transformation p = p = v + /v,v e £(Z), < \2.35\) and < \2.50\) 
change into 

s-RLV n (v, q, r) = s-RLV n (v + jv, q, r) = 0, (2.70) 
RLV n {v,q,r) = RU/ n {v + /v,q,r) = 0. (2.71) 

3 Algebro-geometric Solutions of Stationary Rela- 
tivistic Lotka-Volterra Hierarchy 

In this section, we present a detailed study of the stationary Toda hierarchy. Our 
principle tools are derived from the polynomial recursion formalism introduced 
in section 2 and a fundamental meromorphic function <p on K n . With the help of 
4> we study the Baker- Akhiezer vector ^, the common eigenfunction of and 
V n (£), trace formulas, and theta function representations of 0, p and q. 
Throughout this paper we have the following hypothesis 

Hypothesis 1. Assume 

p — p^ 0,1, p(n),q(n),r(n),p(n),q(n) E £(Z), n e Z, (3.1) 

and n G No in A2.35\) fixed. 

Hypothesis 2. The affine part of the algebraic curve associated with the n-th 
equation in the stationary relativistic Lotka-Volterra hierarchy < \2.35\) . which takes 
the form of 

2n+l 

K n : T n (z,y) = y 2 +4R 2n+2 (z) = y 2 - J] (z-E m ) = 0, {£ m }m=o Q C\{0} 

m=0 

(3.2) 

is nonsingular. That is, 

E m ^E m > for m ^ m ,m,m — 0, 1, 2,- • • , 2n + 1. (3.3) 

lC n defined in (13.21) is compactified by joining two points Poo±, Poo+ 7^ -Poo- at 
infinity, but for notational simplicity the compactification is also denoted by fC n . 
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One can introduce the complex structure on fC n to yield a compact Riemann 
surface of genus n ||2T1 1221 |23l . This is a hyperelliptic Riemann surface. It can 
be regarded as a double covering of sphere surface Coo = C |J{oo}. Points P on 
fcp\{Poo±} are represented as pairs P = (z,y), where y(-) is the meromorphic 
function on JC n satisfying TjyZ, y) = 0. 

We write 

n n 

B n {z) = (-q + ) U(z-H), C n (z) = (-f) lliz-Uj), v, e £(Z). (3.4) 

3=1 3=1 

Next we lift the point fiji v j from the Coo to the compact Riemann surface K, n , 

h = (Pji -^A n+1 (/j,j, n)), Uj = (uj, 2A n+1 (uj, n)), j = 1, 2,- • • , n (3.5) 
and introduce the points P ,± on K-n 

2n+l 

4 



j 2n+l 

P ,± = (0,±2i n+1 (0,n)), A 2 n+l (0,n) = - ] [ E m . (3.6) 



m=0 



We introduce the holomorphic sheet exchange map on K n 

*: K n ^K n P = (z,y)^P*=(z,-y) P^ ^ P^ ± = P^. 
Next one can define a fundamental meromorphic function <p on /C„ 



|y- A n+l {z,n) zC n (z,n) 



The divisor (0(-)) of is 

n)) = V Po + o [n) - £> Poo+ £( n) , (3.8) 

where we abbreviated 

= (/ii(n), /i 2 (n),- • • , jU n (ra)) , z>(n) = (vi(n), u 2 (n),- ■ ■ , u n (n)) . (3.9) 

Given </>(•, n), the stationary Baker-Akhiezer vector n, ra ) on is men de- 
fined by 

* (P '^ n ' no)= U^e,n,n ) J' (3 - 10) 
13 



{n"'=n +i [Pi 71 ') ~ z + <?(™')0 ( P ' n ')) > n > n + 1, 
1, n = no, 

nXn+i ~ 2 + <?(™')0~( p n ')) > n < n - 1, 

(3.11) 



^ 2 (P,e,n,n o ) = r 1 0(^,no)< 



1, 

n no 



r(n )z 
W+l \4>-(P,n) 



r(n )z 



1 



-1 



^W=n+1 \4>-(P,n) 

Some properties of 0, are discussed in the following lemma. 



n > n + 1, 

n = no, 

n < no — 1. 

(3.12) 



Lemma 1. Suppose p, q satisfy the n-th relativistic Lotka-Volterra hierarchy ( 12.351 ) 
and A3.1\) - A3.l2\) holds. Let P = (z, y) G /C„\{P 00± , Po,±}, (n, n ) G Z 2 . TTien 
satisfies the Riccati-type equation 



as well as 



(p - z)4> + 



0(P)0(P* 



0(P) + 0(P*) = 

4>(P) - 4>(P*) 

Moreover, the vector \I/ satisfy 



' = rz + 4>~ 

-zC n {z) 
B n {z) ' 

— 2A„ + i 

y 



BJz) 



U(£)*-(P) = \P(P), 
K(0^(P) = (1/%*(P), 

tf 2 (P, £, n, n ) = r'0(^ n)*i(P, f , n, n ) 



(3.13) 

(3.14) 

(3.15) 
(3.16) 

(3.17) 

(3.18) 
(3.19) 
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Proof. Using (12321) (12341) (12361) (12371) (12391) (13771) . we directly calculate 



(p — z)<f) + q<j)<f) — fz — <p 

= l/{B n B-)[B-{p - z){-y - A n+1 ) + q(-y - A n+l ){\ - A; i+l ) 

1 11 (32 0) 

- fB n B~ - {-y - A n+1 )B n ] 
= 0. 

Equalities (13.14l) - (13.16l) hold from (13.71) . Next we use induction to prove (13.191) . 
Obviously, for n = n , we have 

tt 2 (P, £, n , n ) = rV(P n )^i(P, £, n , n ) (3.21) 

from (13.1 II) and (13.121) . In the case k > n , assume (13.191 ) holds for n = no, n + 
1,- ■ • , k - 1. From the definition of \& 2 in (13.111) (13.121) . one finds 

f tt 2 (P, f , fc, no)/*i(P, £, fc, %) = ^x [tf 2 (P, f , fc - 1, no)/*i(P, £, fc - 1, n )] 

^)/(,~-, + #-(P,,))] 

= £ x r x 0(P fc - 1) x [(^^y) / (p - * + <?<H^ fc ))] 
= (fz + 1) / (p - z + #"(P, k)j . 

(3.22) 

The Riccati-type equation (13.131) shows that 0(P, fc) satisfies 

4>(P, k) = {fz + l)/(p-z + #"(P, fc)) . (3.23) 

Noted that ^^(P, £, k, n )/^i(P, £, k, n ) and 4>(P,k) take the same value at 
n = n and (13.221 ) (13.231 ). we have (13.191 ) for n = k. The proof of the case n < n 
is similar with n > n . From (13.111 ) (13.121) . one finds 



(3.24) 



2 
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and 



vjr 2 = + 1 



where we used (13.191 ) in the last equalities of (13.241 ) and (I3.25I ). (13.181 ) comes from 
0771 ) and 0.19ft by direct calculation. 

□ 



Combining the polynomial recursive relation defined in section 1 with 0.41) 
then yields the following trace formula for ag, and bi in terms of symmetric func- 
tions of fj,j and Uj, respectively. 

Lemma 2. Suppose p, q satisfy the n-th stationary relativistic Lotka-Volterra sys- 
tem ( 12.351) and ( I3.il) . Then we have the following trace formula 

n 

p + — q + f — q ++ f + — q ++ /q + + 8\ = — [Aj, (3.26) 

5=1 

p — g + f — r~/r — qr~ + Si — — Vj. (3.27) 

5=1 

Proof. Comparing the coefficients of z n ~ l of B n (z),C n (z) in (12.601) (12.611) and 
0741) the yield the above formulas 07261) 07271) . 

□ 

Next we turn to study the asymptotic behavior of </>, \&2 in the neigbor- 
hood of Poo± and P ,±. This is a crucial step to construct the algebro-geometric 
solutions of relativistic Lotka-Volterra hierarchy. 

Lemma 3. Suppose p, q satisfy the n-th stationary relativistic Lotka-Volterra sys- 
tem ( 12.351) and 0771) . Moreover, let P = (z,y) G JC n \{Poo±, Po,±}, (n, n ) G 
Z x Z. TTien defined in ( 13.71) /ias the following asymptotic property 

j( P)= j(? + )~ 1 r 1 +(((? + /g-i)^)/?) + +o(c) « p^Poc+, (328) 



-f+(pf -p^)C + C(C 2 ) as -P -> -Poo-, 
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where we use the local coordinate z = ( 1 near the points Poo± 

4>{P)-- 



Cn/{Y[ 2 Z=lE m )C + 0{e) as P^P ,+, 



(m=oEm)/b n + 0(Q as P^P 



(3.29) 



where we use the local coordinate z = ( near the points Po,±. 

The components ^i, ^2 of the Baker-Akhiezer ^ have the following asymp- 
totic properties 



(q+(n)v + (n))/{q+{n )v+{n )) + O{() as P -> P^, 



(_l)n-no x C no_n (1 + 0(0) 



a* P ->■ Poo_, 
(3.30) 



*f P t , fv+(n)/v+(n ) + 0(0 as P -+ P , + , 

Wi(P, t, n, nn) = < rr 2 "+ 1 p (3.31) 

\np-q Ik ^)(n,n ) + O(0 as P -+ P ,_, 

^ rp , 1 = r i (i + w/ (rwwr 1 (i + o(0) <w p^Poo + , 

21 ' ° J 5 \(-l)n+i-n x f(n)C°- n (1 + 0(0) as P^Poo-, 



y 2 (P,Z,n,n )=C 1 X< 



{ [(v + (n)c n ) I (n^o 1 ^ + K))]C + 0(C 

T-r2n + l 7-1 



(3.32) 



as P — > Po,+, 

))/&n + 0(0 



where v are defined in Remark\3\and 

Un'=n + 1 f( n ') U > n 0, 

r(/)(n,no) = <l n = n , V/ e £(Z) 



as P P 0) . 
(3.33)' 



(3.34) 



Proof. The existence of the asymptotic expansion of <fi in terms of the local coor- 
dinate z = C 1 near Poo±> respectively, ( = z near P 0> ± is clear from the explicit 
form of in (13.71 ) (I3.8I ). Assume has the following asymptotic expansions 



^C^ + ^o + ^iC + O^ 2 ), as P^P C 



oo + ; 



0o + <M + <M 2 + O(C 3 ), 



as P -»■ P. 



(3.35) 
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Inserting the asymptotic expansions (13.351 ) into the Riccati-type equation (13.131) 
and comparing coefficients of powers of £, which determines the coefficients fa 
in (13.251) . one concludes (13.281) . Insertion of the polynomials B n , C n defined in 
(12.601) (12.611) into ((3771) then yields the explicit coefficients in (13 29b . Next we 
compute the asymptotic expansion of Noted the definition of in (13.1 II) . we 
first investigate the expression p — z + qfa . With the help of (13.281) (13.291 ). one 
finds 



p — z 



'{q + v+)/{qv) + 0{Q, as 

-C _1 + 0(1), as 

v + /v + 0(()i as 

p-q U ™=l Em +0((), as 



P P<x+i 
P -)• Poo-, 



(3.36) 



which give rise to (13.301) (13.311) . Obviously, the exact n poles of p — z + g</>~ 
in K, n \{P oc ±, Po,±} coincide with the ones of fa and there are n + 1 poles of 
p — z + g</>~ in K, n . However, it is easy to know that ^ is a meromorphic function 
on fC n from the definition (13.1 II) and the meromorphic property of <p. Therefore 
the meromorphic function p — z + qfa possess exact n + 1 zero points on K, n . So 
we need some deformation of function p — z + qfa. Using (12.641) (12.671) (13.21) and 
(13.71) . one may calculate 



p — z + 



p — z + q 



B 



n 

(p - z)B- + q(\y - A' 



qA n+1 + B n + \qy 

B n , + 



+ 9 



I7/ 2 - A 2 



^ B- [\y - A n+1 



B r 



zC n 



B 



kv-A 



n+1 



(3.37) 
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Thus, 

p-z + #"(P) m |^0(1), (3.38) 

which shows /i,-, ^ = 1,- • ■ , n are n zero points of function p — z + q(f>~(P). 
The remaining one zero point is denoted by P$. Finally, (13.321) and (13.331) follows 
from (13191) and (|3T281) - (|33TT) . 

□ 

We choose a fixed base point Q on /C p \{Po )+j ±, Poo±}- Let u}p^ Poo be a 
normal differential of the third kind holomorphic on IC p \{P oc+ , Po,+} with simple 
poles at Poo + and P and residues -1 and 1, respectively, that is, 



(3) C^Q J (-C- 1 + 0(1)K P->Poo+ „,„ 

Up °< +P ~ + 1 (C" 1 + 0(i))dC P ^o,+ 



and 

,(3) 

y 



Wp , p 

00 + -* oo - 



i=i 

be a normal differential of the third kind holomorphic on /C P \{P 00+ , Poo-} with 
simple poles at P^ and P 0O _ and residues 1 and -1, where the local coordinates 
z = C -1 for P near P^ , z = ( for P near P , and {Aj}j=i r . jP , {X-}j=i r . jP are 
constants uniquely determined by normalized process. 
Moreover, 

'"2U. + = °» J = 1"-,P, (3-41) 
(3) C^o / \,f eo._ \ , n , P P 





p 



Uc + u; +o(c) p-iic- (3 - 42) 



/ = -lnC + ^o + 0(C) P^P , (3.43) 

where we choose a homology basis {dj,bj}^ =1 on /C n in such a way that the 
intersection matrix of the cycles satisfies 

ajobk = 6j tk , ajoa k = 0, bjob k = 0, j,k = l,---,n (3.44) 
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and e ,±, d £ C. One easily verifies that dz/y is a differential on /C n with zeros 
of order p — 1 at Poo± and hence 



-d*, j — 1, 



form a basis for the space of holomorphic differentials on /C n . Introducing the 
following invertible matrix Cj^ £ C 



C — {Cj t k)j,k=i,-,pi Cj,k — I Vji 

J a k 



(3.45) 



c(fc) = ( Cl (fc),--- )Cp (fc)), c,(fc) = (C- 1 ),^ fc = l,---,n. (3.46) 

It's easy to show that the normalized holomorphic differentials {uj}j=i,-, P can be 
written into 



22 c j(l)vi, / Vj = 8j,k, j,k = l,---,n, 



(3.47) 



u = (wi,- ••,£*;„). (3.48) 
Assume 77 £ C and |^| <min{|_E' |~ 1 , |-Ea| _1 , l-E^I -1 ,- • ■ , l-EWil -1 } and abbrevi- 



ate 



Then 



where 



Similarly, 



where 



El = (Eq, Ei,- ■■ ,E : 



2n+l; 



f 2n+l 



V2 +oc 

l(l-E mV )\ =Y,^E)v k 



k m=0 



fc=0 
2n+l 



c (E) = 1, Cl (£) = - ^ £ m , etc. 



m=0 



'2n+l 



V2 +co 
1(1-^)1 = J>Clfo* 



. m=0 



fc=0 
2n+l 



co(^) = l, Cl (£) = -^£ 



etc. 



m=0 
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Obviously, 

+00 

y(P) = TC n - 1 Y. c ^ k 

(3.49) 

as P -»■ Poo ± , 2 = C" 1 



In the following it will be convenient to introduce the abbreviations 

z(P,Q) = E Q -A Qo (P)+a Qo (V g ), Pg/C p , Q = {Qt,- ■ ■ , Q p } E Sym n (/C p ) 

(3.50) 

where Sq o is the vector of Riemann constants and the Abel maps Aq q (-), Oq (-) 
are defined by (period lattice L p = {z E 7L 9 \z_ = n + mr, n, m E Z 9 }) 

P^A Qo (P) = (A Qihl (P),..,A QotP (P))= ( r / wp) (3 ' 51) 
and 

a Qo : Di W (X^) J(/C P ),P H- a Qo (P) = £ P(P)^ (P). 

pg/c 9 

Theorem 4. Suppose that p, q satisfy the n-th stationary RLV hierarchy, let P E 
fc n \{Poo±, Po,±} ( n ; n o) £ is non-special. Moreover, 



%(P,A(n ))) V" '" U/ ^Qo 



*i(P, n, n ) = C(n, n o)-^j-^hj—^ x exp ( (n - ra ) / u P|( „, )Pso _ 



(3.53) 

9 2 (P, n, n ) = r 1 x C(n)C(n, n ) §^$^% 

9{z{P, fj,{n Q ))) 

/ nP pP \ 

x exp I / w£ } P + (n - n ) / wp„( n )P K 
flfetPoc-./Kno))) 



C(n,n ) 



^U(Poo-,/i(n))) 
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and finally p, q are the form of 

v += l - (-A 3 - A+ - 6 1 + 1 - Ax ± ((A 3 + A+ + - 1 + A x ) 2 + 4A 2 ) 

(3.55) 

g+ = A 2 /p + + 1 

= 2A 2 (-A 3 - A+ - 6i + 1 - A x ± ((A 3 + A+ + 5 1 - 1 + A x ) 2 + 4A 2 ) 
+ 1. 



(3.56) 



Here 



/ ( 9{z{Poo+i M n )) + <±0\ 

i = V A - - > Cjikjd^ In — — a=0 , (3.57) 

A 2 = Koo+ - > ^(n)^ In — — — r w=0 , (3.58) 

^ v^U( p oo+,M^)) + w)y - 



.7 = 

gU(Poo-,£(7i))) flUCPoo-HftM)) a 2 
3 0U(Poo-,Mw))) 0U(^oo+^W)) 5i 

and 5i, a 2 , {Aj}j=i r . in G C m rt:?.:?2D . 



(3.59) 



Proof. The proof that the divisor T>^ n ) is non-special see Lemma [8l where t r 
is regarded as a parameter. Hence the theta functions defined in this lemma are 
meaningful and not identical to zero. Obviously, 

is holomorphic function on compact Riemann surface /C„(Riemann-Roch Theo- 
rem [14]). So it is a constant C(n) related to n and 0(P, n) has the form (13.521) . 
One have the following expansion as P — >■ Poo + , (z = 



9(z(P,u(n))) _ ejzjP^ujn))) 
9(z(P,Hn))) BiziPoo+^in))) 



x 



j'=i 



5 



g(g(Poo+, j>M) +^) 

(9wj ^(z(P 0O+ ,/i(n)) + w) 



In 



(3.60) 



|.=oC + o(C 2 ) 



22 



Then as P — > Poo + , 

- ai ° W ^U(P 00+ ,Mn))) 

1 ~ l^ c i n hr~ ln wn~5 — \ \ , — r U=oC + o C 



X 



XT 1 (l + wC + O^ 2 )) 



where a±, Koo + G C are constants generated in the limit procedure. In another 
way, meromorphic function has the asymptotic expansion (13.281) . Comparing 
the coefficients of (13.281) with the ones of (13.611) then yields 



(q] = " c,,, Swl' (3 ' 62) 



and 

- l) p + 

^ 3 / g^^^n))-^) ^ . (3.63) 
- Koo+ - ^ c ,(n)_ In ^ U(Poo+ ^ W)+ ^J l^o, 

where u = (ux, oj 2 , - ■ ■ , w n ). Similarly, as P — » Poo-, one finally derives 

- ~ ^ ^U(Poo-,^))) „™ 

- r = a 2 C (n — - — t^tt, (3.64) 

where a 2 G C. Noted the assumption (12.11) (13.11) that q(n) ^ 0, f (n) 7^ for all 
n G Z, we conclude 

-+-_ a 2 gU(Poo-^H))^U(Poo + ^H)) 

^ S^^Poo.,^)))^^^,^))) ( ' ^ 

from (13.621) and (13.641) . Let us consider the trace formula (13.261) . After a standard 

residue calculation at Poo± lfl8l [T9l . £)?=i/4j has the following theta function 
representation 

^ = J - L ln ^ U (P 00 _,^ ))+ ^ ) J 1^ (3 - 66) 
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Then p, q satisfy the following equations 



(q ++ /q + - 1) p + = A 



2) 



(3.67) 



- q + r = A 3 , 



(3.68) 



Plugging (13371) (I3T68T) into (f3T26l) then yields 

(p+) 2 + (A 3 + A+ + Ai - 5 1 + 1) p+ - A 2 = 0. 



(3.69) 



The two solutions of (13391) are (13351) . (13361) follows from (13371) and (13351) . 

□ 

Remark 5. 77ze theta function representations of^i, ^ ZM ( 13.531) ( 13.541) z's merely 
a formal expression and they are of little use in the process of solving the station- 
ary relativistic Lotka-Volterra system ( 12.351) . In fact the expression ( 13.531) itself is 
indeterminate as Pjj is unknown and may relate to the lattice variant n. 

4 Algebro-geometric Solutions of Time-dependent Rel- 
ativistic Lotka-Volterra Hierarchy 

In this section, we mainly derive the algebro-geometric solutions of time-dependent 
relativistic Lotka-Volterra hierarchy defined in section 2 by extending the method 
employed in the section 3 to time-dependent cases. 

Hypothesis 3. Assume 



and (n, t) E Z x R. 

Throughout this section we suppose Hypothesis 2 and Hypothesis 3 holds. 
The basic algebro-geometric initial value problem is that if we consider a so- 
lution p l {n),q l {n) of the n-th stationary relativistic Lotka-Volterra system s- 



=0, associated with the hyperelliptic curve K n and a corresponding of the summa- 
tion {<^}" =0 C C, then we construct a solution^, q of the r-th time-dependent rela- 
tivistic Lotka-Volterra flow RLV r (p, q) — satisfying p(n, t 0>r ) = p 1 (n), q(n, t 0>r ) 

= g x (n) for some t , r E R and any n E Z. We shall use the notation V r , A r+ i, B r , 



p = p^0,l p(;t),q(;t),p(;t),q(;t)f(;t)E£(Z) 

p(n, •), q(n, -),p(n, •), q(n, •), f(n, •) E C^R) 



(4.1) 



RLV p {p\q x ) 
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C r , D r+ i, at, bi, &i, de, 8 \ in the r-th time-dependent flow to distinguish A n+ i, B n , 
C n , D n +i, a>£, be, q, di, Si in the n-th stationary relativistic Lotka-Volterra system. 
The algebro-geometric initial value problem discussed above can be summed up 
in the form of zero-curvature equation 

UtriM + U(£,tr)V r {Z,t r ) - V r + (Ur)U(Ur) = 0, (4.2) 

U(Z, t ,r)V n (£, V) - t 0ir )U{£, to,,) = 0, (4.3) 

Considering the isospectral property of the Lax operator L corresponding to U, 
we may impose more strong condition on equation (14.31) . which is 

U(Z,t r )V n (Z,t r ) - V+(£,t r )U(£,t r ) = 0, tr G R. (4.4) 

For further reference, we recall the relevant quantities here: 

m = ( p z z ? ) , (4.5) 



where 



v ^ = { -Cw ) ' <4 ' 6) 
mo J few f:« V (4.7) 

n+1 

A n+1 (z) = ^d n+l ^z l , (4.8) 

1=0 
n 

B n (z) = Y,bn-ez i , (4.9) 

1=0 



C n (z) = Y,c n -iZ l , (4.10) 

£=0 
n+1 

D n+1 (z) = ^2d n+1 _ £ z e , (4.11) 



£=0 
r+1 



ir+l(«) = ^ Or+l-^, (4.12) 



e=o 
25 



B r {z) = Y}r-iz\ (4.13) 

r 

6 r (z) = Y,c r -iz £ , (4.14) 

1=0 
r+1 

D r+1 (z) = J^2 r+ i_^, (4.15) 

Here {a^g, W^o and {^}|=o> {%=o are 
defined by (I2.2I) - (I2.6I) corresponding to different constants 8i and respectively. 
Explicitly, equation (14.21) and (14.41) are equivalent to 



- Pu = {p~ z)A r+l + zqC r - (p - z)A r+ i - zrB r , 


(4.16) 


~ Qt r = (p - ^)5 r " - - gi r+ i - B r , 


(4.17) 


- f tr = fA~ +1 + C~ - (p - z)C r + fD r+1 , 


(4.18) 


= zfB~ - D~ +1 - zqC r + D r+1 , 


(4.19) 


= (p - + zgC*" - (p - z)i n+1 - zfB n+1 , 


(4.20) 


0=(p-z)B~- qDn+i - A+i - B n , 


(4.21) 


= ri~ +1 + C~ -{p-z)C n + rD n+1 , 


(4.22) 


= 2X5" - D~ +1 - ZqC n + D n +1, 


(4.23) 



respectively. In particular, (12.671) holds in the present t r -dependence setting, that 
is, 

R 2n +2{z,t r ) = -(A n+1 (z,t r )) 2 - zB n (z,t r )C n (z,t r ). (4.24) 

Obviously the algebraic curve defined in (14.241 ) is n -independent and may depend 
on the parameter t r . In fact we can prove 

d tr R 2n +2{z,t r ) = 
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under the initial value condition ( 14.41) , which means R 2n +2(z, U) is t r -independent 
(see lemma [5]). We write 



B n (z, n, t r ) = (-q + (n, t r )) JJ (z - fj,j(n, t r )) , 



3=1 



(4.25) 



C n (z,n,t r ) = (-f(n,t r )) (z - Uj(n,t r )) , 

3=1 

and jjLj{n,t r ), Uj(n,t r ) G ^(Z), (n, t r ) eZxR. 
As in the stationary context (13.51) (13.71) we introduce 

p,j(n,t r ) = (fij(n,t r ), -2A n+1 (iij(n,t r ),n,t r )), 

Pj(n,t r ) = {v j (n i t r ) i 2A 7H . 1 (j/j(n,t r ),n,t r )), (4.26) 

j = l,2,---,n, (n,tr) GZxR. 



on /C n and define the following meromorphic function 0(-, n, t r ) on /C 

- A n+1 (z,n,t r ) zC n (z, n, i r ) 
B n (z,n,t r ) \y + A n+1 {z,n,t 



<j){P, n, t r ) = — — = 5 — , P = (z, ?/) e £ n . 



(4.27) 

The divisor (</>(•)) of is 

(0(-, n)) = V Po + o(n,t r ) ~ ^F 0o+ A(n ) tr)> ( 4 - 28 ) 

and the time-dependent Baker- Akhiezer vector is then defined in term of by 



*i(P,^,n,n ,tr.,to,r) = exp ^ (a. + i(z, n , s) + P r (z, n , s)(j)(P, n , s) ) ds 

p(n ,t r ) - z + q(n')<f>~(P,n ,t r )\ , n>n + l, 

n = no, 

JT^ =n+l (p{n\t r ) - z + q(n\t r )^-(P,n\t r ) s j , n < n - 1, 




(4.30) 
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$> 2 (P,£,n,n ,t r ,t 0)r ) = exp / ( A r+ x(z, n , s) + B r (z,n , s)0(P,n o , s)) ds 



n: 



x f 1 x 4>(P, n , t r ) < 



n'=n +l ^-(P,n',t r ) 



n > n + 1, 
n = n , 



(4.31) 



n 



"0 



-1 



/^n =n+l \<^-(P,n ,t r ) 

P = (z, y) G /C n \{P ,±, Px,±}, (n, t r ) G Z x R. 
One observes that 

^i(P, n, n , t r , t , r ) = ^i(P, n , n , t r , t 0>r )^i{P, n, n , t r , t ,r) 
P = (z, y) G /C n \{P ,±, Poo±}, (n, n , * r , V) e Z ' x ^ 



(4.32) 



(4.33) 



The following lemma shows the properties of 0, and \l/ 2 as discussed in the 
stationary case. 

Lemma 4. Assume Hypothesis 2 and Hypothesis 3 hold and suppose p(n, t r ), q(n, t r 
satisfy ( 14.21) H4.4\) . In addition, let P = (z, y) G K, n \{P OQ ±, Po,±}, (n, n , t r , t 0jr ) G 
Z 2 x I 2 . Then 0(P, t r ) satisfies the following equations 



as well as 



(p - z)(j)(P) + #(P)0~(P) =rz + 0~(P), 

V(P) = " 4+10(P) " ^W(P) - Pv0 2 (P), 



0(P)0(P*) = ■ 

0(P) + 0(P*) = 
0(P)-0(P*) 



y 



BJz) 



The vector \l> satisfy 



u(t)*-(P) = *(P), 
K(0^(P) = (1/2)^(P), 
2 (P, £, n, n ) = r'kP, n^P, £, n, n ), 



(4.34) 
(4.35) 

(4.36) 

(4.37) 
(4.38) 

(4.39) 
(4.40) 
(4.41) 
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*tr{P)=V r (ZMP)- ( 4 - 42 ) 

Moreover, Moreover, as long as the zeros of fa (no, s) ofB n (-, n , s) are all simple 
and nonzero for all s G £1, Q C R is an open interval, is meromorphic on 

fc n \{Po,±, Poo±} far (n,t r ,t 0t r) eZxfl 2 . 

Proof The proof of (14.341 ), d4.36M4.411 ) is the same with lemma [B where t r is 
regarded as a parameter. From ( |4.34b we have 

Pt r 4> + (p - 2)0 + % r 00" + #< r 0~ + = h r z + (4.43) 

that is, 

\p-z + q<j>~ + (q<f> - l)^ - ]^ = f tr z - p tr (f) - q tr 4><fi~ '• (4.44) 
Using d4.16M4.18ft . one finds 

[p - z + qfa + (q<p - l)S~] (zC r - D r +i4> ~ A r+1 (j) - B r 4> 2 ^ 

= + (# - l)S~] (zC r - A.+10 - - 5 r 2 ) 

= (p - 2 + qcj)~)zC r - (fz + <fi~)D r+1 - (fz + 0~)A r+ i - (fz + <fi~)B r <f)~ 

+ (# - i)zd r - - (# - i)£>- +1 0- - (# - i)A; +1 fa - (# - 1)40 2 

= (p - z + q(j)~)zC r - (fz + 0-)-D r+ i - (fz + 0~)A r+ i - (fz + <fi~)B r <f)~ 

+ (# - i)*d r - - (# - - (# - i)i r - +1 r - f«J r -0- + s r (p - z)h> 

= [fA~ +1 - C~ + (p - z)(7 r - f.D r+ i]z + [(p - 2)i r 7 +1 + zgC^ - (p - z)A-+i 
- zrB r+1 ]<j> + [(p - 2)57 - q~D~ +1 - qA r+ i - B r ]4>4>~ 

+ [(A; +1 - A r+1 )(fa - (p - z)<j> - q<$T + fz)} + [zfB- - D; +1 - zqC r + D r+ 
= f tr z - p tr 4> - q tr H~ ■ 

(4.45) 

Considering (14.44b . one may conclude 

\p-z + qfa + (# - l)S-\ (fa - (zC r - D r +i4> - - 5 r 2 )) 

= 0. 
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Thus, the expression 

{YZ' =no+ xB{z,n ,t r ), n>n + l, 
1, n = n , 

(4.47) 

where 

x 1 — q(n,t r )<j)(z,n,t r ) , _ _ A 

B(z, n, t r ) = qK '_ rm - -, (n, * r ) GZxi (4.48) 

p{n,t r ) — z + q{n,t r )<p {z,n,t r ) 

Obviously, the expression <p tr — (zC r — D r +i4> ~ A r+ i4> — B r cj) 2 ) in the left-hand 
of (14.471) is a meromorphic function on Riemann Surface K, n and its order is finite 
at Poo-- However, the order of the right-hand in (14.471) is O(z n ~ no ) as P near the 
point Poo_. Hence, taking n sufficiently large, then yields a contradiction at the 
both bides of (14.471) unless C = 0. Therefore we have (14.351) . From the definition 
of ^i(P, £, n, n , t r , t 0t r) in (|4.30l) it is easy to check 

(p - z + q4>~) tr = Pt r + qtr<jr + # t ~ 
= — (p — z)A r 7 +1 - zg(7~ + (p — z)A r+ i + zfP r 

X (~(P - ^)^r~ + + rfr+1 + #r) 0" (4 49) 

= (p - z + q<jr\ [A r+1 + B r 4> - A~ +1 - B~(j> 

If we note 



T(P,n ,t r ,V) =exp (/ (A+i(^"o,s) + P^(z,rio,s)0(P,n o ,s)) rfsj , 



(4.50) 
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then 

*i,t, = ( T f[ (p-z + q4>-)(n) 

\ n'=n +l 
n 

= T *r n (p- z +v<f>~)( n ') + t (^ _ ^+#~)*r( n ') n 

n'=no+l n'=?io+l n" 

n n 

= T JJ (p-z + ^-)(n')[i r+ i(no,t r )+ J] (i r+1 (n) + 

n =no+l n =no+l 

i r - +1 (n')- B r " (n))] 

= (Ar+l + S r )*l = A-+l*l + f#r*2, 

as n > no + 1. 

(4.51) 

The proof for the case n < n — 1 is similarly with (14.511) . Using (14.341) and (14.511) 
we have 

* 2A . = r 1 (0^*i + 0*i,*j 

= r 1 - A-+10 - - 5 r 2 + (A r+1 + B r+ i0)0) * x (4.52) 

= r\zd r - 4+0*1 = £4*1 - 4 + i*2. 

To illustrate that (•, £, n, n , t r , t ,r) is a meromorphic function on /C n \{P ,±, P>o±} 
we only need to prove T is a meromorphic function on /C re \{Po,±, Poo±}- Taking 
into account (14.261) ( 14.271) (14.541 ), one may derive 

B r (z, no, s)<j)(P, no, s) ~ d s In (S r (^, n , s)) + 0(1) (4.53) 

as P — > Hj(n , s). Hence we conclude T is a meromorphic function on /C n \{P ,± ; -Poo±} 
with the help of analysis the asymptotic behavior of the possible poles {fij(no, s)} 

of the function A r+ i(z, n , s) + B r (z, n , s)<j)(P, no, s). 
Next we consider the t r -dependence of A n+1 , B n , C n . 

Lemma 5. Assume Hypothesis 2 and Hypothesis 3 hold and suppose p(n, t r ), q{n, t r ) 
satisfy (\4.2\) (\4.4\) . In addition, let (z, n,f r )eCxZxK. Then, 

B n>tr = -2A n+1 § r + (i r+1 + D r+1 )B n , (4.54) 



) 
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A n+htr = z(B r C n - B n C r ), (4.55) 

C nA = 2z~ 1 A n+1 d r - C n (A r+1 + D r+1 ). (4.56) 
In particular, (\4.54\) - ft4.56\) is equivalent to 

V n ,t r = [V r ,V n ]. (4.57) 
and the algebraic curve defined in A4.24\) is t r -independent. 
Proof. To prove ( 14.541) one first differentiates equation (14.381) 

kXP)~kXn = Z] &^- (4-58) 

n 

The time derivative of given in (14.341) and (14.371) (14.381) yields 

k(P)-k(P*) 

= § r (4>(P) + 0(P*))(0(P) - 0(P*)) + (A r+1 + b r+1 )(4>(P) - 0(P*)) 

- ^ [A-r+l + L> r +l)-=- 

B n R n 

(4.59) 

and hence 

B njtr = -2A n+1 B r + (X+i + B r+1 )B n . (4.60) 
Similarly, starting from (14.371 ) 

<f> tr (P) + 4> tr (P*) = 2zC r + (A + D r+1 ) Z -^ - B r (^ + ^P) (4-61) 

B„ B„ D r , 



II 



yields (14.551) . Differentiating the equation ( 14.231) then yields 

z h r B~ - zrB~ tr - D~ +l tr - zq tr C n - zqC nttr + D n+lttr = 0. (4.62) 
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Using (TX221) (14T71) (|4~T8T) <l4~T9l) (|4~2TT) (14321) d423l) (14341) and (14351) we have 

= -zB~[fA; +1 + d r 7 - (p - z)<% + fb r+1 ] - zf[-2i; +1J §- 

+ {A~ +l + D~ +1 )B~] + z[C n B r — B n C r — C n B~ + B~C~\ 

+ - z)B~ - qE>r+i ~ <Mr+i - B r ] - zqC nA 

= -zB~[-(p - z)C r + fb r+1 ] + zr[-2A- +1 B- + D; +1 B~] 

+ z[-B n C r - C~B~] + z(7 n [(p - z)B~ - qD~ +l - qA r+1 ] - zqC n , tr 

= z[(p- z)B~ - B n ]C r + z[-2rA~ +1 - C~ + {p- z)C n }B; 

+ z[fB~(D~ +1 - D r+ i) - qC n D~ +1 ] - zqC n A r+1 - zqC n±r 

^ (4.63) 

= ZqC r [2A n+1 + + A»+l)] + (^n+l - 

+ [(— -D„+i + D~ +1 + zqC n ){D; +l - D r+l ) - qC n D; +l ] 
— zq~C n A r+ \ — zqC n> t r 

= 2zqA n+1 C r + (I re+ i - I- +1 )(^rS- - zg(7 r ) + [(-D r +i + 
D~ +1 + zqC n )(D~ +1 - D n+ i) - qC n D~ +1 ] - zqC n A r+ i - zqC nA 

= 2zqA n+1 C r + (A n+1 - A- +1 )(D; +1 - D r+1 ) + [(-D n+1 + 
D~ +1 + zqC n )(D~ +1 - D r+ x) - qC n D~ +1 ] - zqC n A r+1 - qC n>tr 

= 2zq~A n+ iC r — zqC n D T+ i — zqC n A r+ i — zqC n ± r = 0, 

which is equivalent to (14.561) . Finally one can directly differentiate (14.241) 

R tr = -2A n+1 A n+Ur - zB nA C n - zB n C n ± r (4.64) 

and insert (14.54b - d4.56b into (14.641) to yield R tr =0. □ 

□ 

Next we derive Dubrovin-type equations, that is, first-order ordinary differen- 
tial equations, which govern the dynamics of jij and Uj with respect to variations 
of t r . We recall that the affine part of K, n is nonsingular if 

{£ m } m=0 ,...,2n+i C C, E m ^ E m > for m ^ m' , m, m = 0, . . . , 2n + 1. 

(4.65) 
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Lemma 6. Assume Hypothesis 2 and Hypothesis 3 and suppose ( 14. 2 D ( 14.41) /zo/rf 
onZx with f2 M CI an open interval. In addition, assume that the zeros Hj of 
B n (-) remain distinct and nonzero onZ x TTien {Aj}i=i,- :« defined in ( 14.2(51) 
satisfy the following first-order system of differential system on Z x 

with 

fij(n, •) e C°°(^, JC n ), j = l,...,n,neZ. 

For the zeros Uj ofC n , identical statements hold with fij and replaced by Uj 
and Q u , etc. In particular, {pj}j=\ r .,ni defined in ( 14.2(51) . satisfies the following 
first-order system on Z x Q v , 

2f- x y(j> j )CJv j ) 

Vj,tr= rrn ~ ~ ( — y 3 = V ' " , ™, (4.67) 

wzY/z 

i>j(n, •) G C 00 ^,^), j = l,...,n,nGZ. 

Proof. It suffices to consider (14.661) for fJ,j,t r - Using the product representation for 
B n in (14251) 

n 

S„ = (-?*■) II(A - A*i) 

5=1 

and employing the (14.261) and (14.541) . one computes 



proving (14.661) . The case of (14.671) for z^ is analogous using the product presenta- 
tion for C n in (14.251) and employing (14.261) and (14.561) . 

□ 

Since the stationary trace formulas for de, b e in terms of symmetric functions 
of Uj and Uj in Lemma [2] extend line by line to the corresponding time-dependent 
setting, we next record their i r -dependent analogs without proof. For simplicity 
we confine ourselves to the simplest ones only. 
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Lemma 7. Assume Hypothesis 2 and Hypothesis 3 hold and suppose p(n, t r ), q(n, t r ) 
satisfy (\4.2\) (\4.4\) . Then, we have the following trace formula 



p+ _ _ q++ r + _ q++ /q+ _|_ ^ = _ ^ ^ (4.68) 

i=i 

n 

p — q + f — r~ /r — qr~ + Si — — Vj. (4.69) 

i=i 

Lemma 8. Suppose p, q satisfy the rt4.il) < \4.2\) and the n-th stationary RLV system 
< \2.35\) . Moreover, let (n,t r ) G Z x K. and Dp,, ft = (/ti,- • • , jl n ) G Sym n {fC n ), 
^vj\L — (Air ' ' 5 v n ) G Sym n (K. n ) be the pole and zero divisors of degree n, 
respectively, associated with p, q and defined according to A4.27\) . that is, 

fij(n,t r )= ^fi j (n,t r ),-2A n+1 (fi j (n,t r ),n,t r )j G K, n , j = l,---,n, 

Uj(n,t r ) = (vj(n,t r ),2A^ 1 (v j (n,t r ),n,t r ) \ G JC n , j = 
Then T>f^ n ,t r ) an dT)o<n,tr) are non-special for all (n,t r ) GZxi 

Proof. We are only to prove the conclusion for T>^ tr y 2}#( n ) is non-special if 
and only if {/ti(n),- ■ ■ ,(x n (n)} contains one pair of {/tj, /t*}. Tlence, £>/t(n) is non- 
special as long as the projection pj of jlj are mutually distinct, Pj(n) =f Hk{ n ) f° r 
j ^ k. If two or more projection coincide for some n G Z, for instance, 

Hh{no) =■ ■ ■ = Vj k (n Q ) = Ho, k>l. 

There are two cases in the following associated with fi . If fi G C\{E ,- ■ ■ , E 2n +i}, 
then A n+ i(fi , n , t r ) ^ 0. It is obvious that /Un(^o> t r ),- ■ ■ , Pj k (n , t r ) all meet in 
the same sheet and hence no special divisor can arise in this manner. If p, equals 
to some E mo and k > 1, one concludes 

i? n (z,n ,t r )=^O((z-£ m0 ) 2 ) 

and 

But one observes R 2 p+2(z,no,t r ) = —A 2 n+1 — zB n C n = O {{z — E mo ) 2 ) . This 
conclusion contradict with the hypothesis that the curve is nonsingular. We have 
k — 1. Therefore no special divisor can arise in this manner. Then we have 
completed the proof. □ 
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Next, Now we turn to asymptotic properties of 4>,^ 1 ,^ 2 defined in (14.271) 
( 14.301 ) and ( 14.311) in a neighborhood of P ,± and P^. 

Lemma 9. Assume Hypothesis 2 and Hypothesis 3 hold and suppose p(n, t r ), q(n, t r ) 
satisfy ft4~2\) fi4~4\). Moreover, let P = (z,y) G /C n \{Poo ± , P }, (n,n ,t r ) G 
Z x Z x R. Then, the meromorphic function on K. n has the following asymptotic 
behavior 



(q+)-i(-i + (((q+/q - l)p)/q) + + 0(() as P P^ + , 
-r + (pr~ — pr)( + 0(( 2 ) as P — > Poo-, 



where we use the local coordinate z = C, 1 near the points Poo±- 

>c n /(Yl 2 :^E m )( + 0(e) as P^P 0l+ , 
^ j " 1 -(nS^)/^0(() as P^P ,_, (4 - 71) 

where we use the local coordinate z = ( near the points Po,±- 

Proof. By the definition of in (I4.27I) the time parameter t r can be viewed as an 
additional but fixed parameter, the asymptotic behavior of remains the same as 
in Lemma [3] □ 

Theorem 6. Assume Hypothesis 2 and Hypothesis 3 hold and suppose p(n, t r ), q(n, t r 
satisfy ft4~2\) fi4~4\) . Moreover, let P G /C n \{Poo±, P 0j ±} and {n,n ) G Z 2 . TTzerc 
tr \ is non-special. Moreover, 



e(z(P,P(n,t r ))) ( [ P , M 
0(z(P,{i(n,t r ))) P Wq 

and p, q are the form of 



0(p, n, t r ) = c(n, t r ) 2^r:{{{ ^ ( / <oU+ ) > ( 4 - 72 ) 



P + (n, t r ) = i (-A 3 - A+ - <5i + 1 - A x ± ((A 3 + A+ + ^ - 1 + Ax) 2 + 4A 2 ) 

(4.73) 

g + (n,t r ) = A 2 /p + + 1 

= 2A 2 (-A 3 - A+ - 5i + 1 - Ai ± ((A 3 + A+ + 5 1 -1 + Ax) 2 + 4A 2 ) 
+ 1. 

(4.74) 



iv 1 
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Here 



A 2 = Koo+ - > Cj n a In — — y a=0 , (4.76) 

Aa = fl(z(Poo-,£(n,t r ))) gU(Ao+, A(w,t r ))) a 2 

9(z(Poo-, jl(n, t r ))) 6'(2(P 00+ , z>(n, t r ))) ai 

and 5i, 02, {Aj} J= i... >n G C in (15.321) . The Abel map evolving with respect to t r 
are 

d d P r 

■z-2Q A' D Kn,tr)) = ^-°LQoA V z(n,tr)) = ~ EE c e (k)6 r _ s c k+s _ p (E) . 

(4.78) 

Proof. The proof of (14.72b - fl4.77b is analogous with theorem HI Here t r can be 
regarded as a parameter. Next we prove (14.78b . 



d d n P 

Q-«Qo,/( X, A(n,tr)) = E / 

°tr °tr ~\ JO, 



i=i 

n 
3=1 

n / n jjf?~^ \ ( - n 

E [j2 C ^ k \(n.L t )) (9 + )" 1 ^r(/ij(^,^))y(/ii(^tr)) JJ (jij ~ Hk) 

j=l \k=l VWK ' T >> ) \ k=l,k^j 

n / n , k— 1 



g(S c ' w ira^j( (r,J ' fe(M ' )) ) 

- E <*(*) E n » ^—7i E i+i- E ^)*^) 

k=l j=l llfe=l,fc^i W y s =0 y=max{0,s-p} 

p r 

- E E c K fc )^-sCfc+ s -p(p). 

fe=l s =0 

(4.79) 
□ 
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5 Appendix: The Lagrange Interpolation Represen- 
tation of — (^ + )~ 1 J B r+ i(/i J (n, t r )) 

We search for the interpolation representation of B r (fj,j(n,t r )) as in the KdV, 
AKNS, Toda cases. Introducing the notation in lfT8l[T9l , 

Sk = {L = (4,- ■ ■ , 4) e N fc |4 <• • ■ < 4 < n}, k = 1,- • ■ n, 
= (-l) fc ]T w,-,^, 

= U = (4,- ■ • , 4) e N fc |4 < ■ ■ ■ < 4 < n, £ m ^ j, m = 1,- • • , k}, 
k = l,- --n-l, j = l,---,n 

and the formula 

k 

^* fc _^)/4 = $]fV), k = 0,---,n, j = l,---,n. (5.1) 

1=0 

Let B n = {q + )~ x B n , 6 S = = 0, 1,- • • , n), then one finds 

n n n 

B n { Z ) = h n~sZ S = Y[ (Z - fXj) = ^P-I(H)^ 

1=0 3=1 1=0 

and 

6 £ = ^(/i), £ = 0,---,n. 
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In the case r < n, 



r r I min{r— s,n} 

B r = ^b r . s z s = ^2\ Yl Cr-s-k{E)b k I z s 



s= 


a 


r 

El 


/ r—s 

P 


s=0 


\k=0 


r 

El 

s=0 


/ r—s 
\k=0 



s=0 \ k=0 



s=0 t=0 

Using (15.11) . we have 



s=0 t=0 



(5.2) 



(5.3) 
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In the case r > n, 

r r I min{r— s,n} \ 

B r (z) = ^2b r _ s z s = E E c r - s -k(E)b k z s 

s=0 s=0 \ k=0 J 

r—n n r r—s 

= E E Cr-s-k(E)^ k (^) Z S + E Cr+ls-k(E)^Mi 

s=0 k=0 s =r-p+l k=0 



r—n n r n 



EE^- 5 -^)^^)^ e Y dr ~ s ~ k ^^ k ^ zS 

s=0 k=0 s =r-n+l k=0 

n r 

^^c r _ s _ fc CE)^ fc (/i)z s 

k=0 s=0 

r n 

^^c r _ s _ fc (E)^ fc (/i)z s 

s=0 k=0 

r n 

s=0 k=0 

r—n / n \ r / n 

E^(£) E^(/^ p ~ fc K~ n ~ s + E c s (e)[Y^m 

s=0 \k=0 J s =r-p+l \ 

r—n r / r—s 

Yc s (E)(B n (\))z r ~ n - s + E c s (E)[YMti) 



^r—s—k 



r-s-k 



s=r— n+1 \k=0 



(5.4) 

Then we have 



r—s 



B r (^)= E d s(E) E^K" 8 "" = E C.CS)^.^- (5-5) 

s=r— n+1 \fe=0 / s=r— n+1 

Combining (15.31 ) with (15.51 ), one finds 

r 

B r (z)= E c a CB)¥?2,Qx). (5.6) 

s=max{0,r— n} 

Hence 

r r I s \ 

Brfa) = E 1-sB^j) = E E ■ ( 5 - 7 ) 

s=0 s=0 \t=max{0,s-n} / 
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